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Abstract. Fully numerical mesh solutions of 2D quantum equations of Schrodinger 
and Hartree-Fock type allow us to work with wavefunctions which possess a very flex- 
ible geometry. This flexibility is especially important for calculations of atoms and 
molecules in strong external fields where neither the external held nor the internal 
interactions can be considered as a perturbation. The applications of the present ap- 
proach include calculations of atoms and diatomic molecules in strong static electric 
and magnetic helds. For the latter we have carried out Hartree-Fock calculations for 
He, Li, C and several other atoms. This yields in particular the first comprehensive in- 
vestigation of the ground state configurations of the Li and C atoms in the whole range 
of magnetic helds (0 < B < fOOOO a.u.) and a study of the ground state electronic 
configurations of all the atoms with 1 < Z < 10 and their ions A + in the high-field 
fully spin-polarised regime. The results in a case of a strong electric field relate to 
single-electron systems including the correct solution of the Schrodinger equation for 
the H^~ ion (energies and decay rates) and the hydrogen atom in strong parallel electric 
and magnetic fields. 



I INTRODUCTION 



Theoretical studies of atoms and molecules in strong external fields are motivated 
by several applications. The latter are e.g. experiments with intense laser beams 
(electromagnetic fields with dominating electric component) and astronomical ob- 
servations of white dwarfs and neutron stars (magnetic fields). The experimental 



availability of extremely strong electric fields in laser beams makes the theoretical 
study of various atomic and molecular species under such conditions very desirable. 
The properties of atomic and molecular systems in strong fields undergo dramatic 
changes in comparison with the field-free case. These changes are associated with 
the strong distortions of the spatial distributions of the electronic density and cor- 
respondingly the geometry of the electronic wavef unctions. This complex geometry 
is difficult for its description by means of traditional sets of basis functions and 
requires more flexible approaches which can, in particular, be provided by multi- 
dimensional mesh finite-difference methods. 

Let us discuss the problem of atoms in a strong magnetic field in more detail. We 
start this consideration with the hydrogen atom which was the first atom whose 
behaviour in strong magnetic fields was investigated (for a list of references see 
[1-4]). In cylindrical coordinates (p, z) its non-relativistic Hamiltonian has the 
form (we use atomic units throughout our work) 

1 / d 2 1 d d 2 m 2 \ m 7 2 2 1 

H = -2\W + ~pd-p + d^~V) +Szl + ^ + ^ P ~r (1) 

where m is the magnetic quantum number, s z is the spin z projection and 
7 = B/Bq, Bq = hc/ea 2 ) = 2.3505-10 5 T is the magnetic field strength in atomic 
units. The magnetic field is parallel to the z axis. The Hamiltonian (1) contains 
two potentials of different spatial symmetries: the spherical-symmetric Coulomb 
term 1/r and the cylindrically symmetric potential of the magnetic field 7 2 p 2 /8. 
When considering the impact of the competing Coulomb and diamagnetic interac- 
tion it is reasonable to distinguish the three different regimes of weak, strong and 
intermediate fields. In the latter case the magnetic and Coulomb forces are compa- 
rable. In the case of relatively weak fields the main features of the geometry of the 
wavefunction are determined by the dominating Coulomb term whereas the effect of 
the magnetic field can be consider as a perturbation of the Coulomb wavef unctions. 
For the opposite situation of very strong magnetic fields and dominating cylindrical 
symmetry the adiabatic approximation [5-7] was the main theoretical tool during 
the last four decades. This approximation separately considers the fast motion of 
the electron across the field and its slow motion in a modified Coulomb potential 
along the field direction. Both early (see [8]) and more recent works [2,9-12] on the 
hydrogen atom have used different approaches for these regimes of the magnetic 
field. All these calculations had problems when considering the hydrogen atom in 
fields of intermediate strength. The detailed calculations of the hydrogen energy 
levels carried out by Rosner et al [2] also retained the separation into different 
regimes of the field strength by decomposing the electronic wave function either in 
terms of spherical (weak to intermediate fields) or cylindrical (intermediate to high 
fields) orbitals. A solution allowing to obtain comprehensive results on low-lying 
energy levels of the hydrogen atom for arbitrary field strengths including the in- 
termediate field regime is provided by the multi-dimensional mesh solution of the 
Schrodinger equation [3]. 



For different electronic degrees of excitation of the atom the intermediate regime 
is met for different absolute values of the field strength. For the ground state this 
regime is roughly given by 7 = 0.2 — 20. For atoms with several electrons there are 
two decisive factors which enrich the possible changes in the electronic structure 
with varying field strength compared to the one-electron system. First, we have a 
third competing interaction which is the electron-electron repulsion and, second, 
the different electrons feel very different Coulomb forces, i.e. possess different one 
particle energies, and consequently the regime of the intermediate field strengths 
appears to be the sum of the intermediate regimes for the separate electrons. 

The fact that most methods have problems in the intermediate field region has 
consequences for the current state of the art of our knowledge on on multi-electron 
atoms in strong magnetic fields. There exist a number of such investigations in 
the literature [6,7,13-22]. The majority of them deals with the adiabatic regime 
in superstrong fields and the early works are mostly Hartree-Fock (HF) type cal- 
culations. There are also several early variational calculations for the low-field 
domain [15,23,24]. HF calculations for arbitrary field strengths have been carried 
out in refs. [2,19] by applying two different sets of basis functions in the high- and 
low-field regimes. As a result of the complicated geometry this approach suffers in 
the intermediate regime from very slow convergence and low accuracy of the cal- 
culated energy eigenvalues. Accurate HF calculations for arbitrary field strengths 
were carried out in refs. [18,20] by the 2D mesh HF method. Investigations on the 
ground state as well as a number of excited states of helium including the corre- 
lation energy have recently been performed via a Quantum Monte Carlo approach 
[22]. Very recently benchmark results with a precision of 10~ 4 — 10~ 6 for the en- 
ergy levels have been obtained for a large number of excited states with different 
symmetries using a configuration interaction approach with an anisotropic Gaus- 
sian basis set [25,26]. Focusing on systems with more than two electrons however 
the number of investigations is very scarce [6,21,27]. In view of the above there is 
a need for further quantum mechanical investigations and for data on atoms with 
more than two electrons in a strong magnetic field. For the carbon atom there exist 
two investigations [6,7] in the adiabatic approximation which give a few values for 
the binding energies in the high field regime and one more relevant recent work by 
Jones et al [21]. On the other hand, our two-dimensional Hartree-Fock approach 
allowed us recently to perform precise and reliable consideration of a series of multi- 
electron atoms for the whole range of the magnetic field strengths from 7 = up 
to ultrastrong fields 7 = 10 3 - 10 4 [18,20,28-31]. 

II TWO-DIMENSIONAL MESH HARTREE-FOCK 

METHOD 

Our calculations for multi-electron atoms in magnetic fields are carried out under 
the assumption of an infinitely heavy nucleus in the (unrestricted) Hartree-Fock 
approximation. The solution is established in the cylindrical coordinate system 



(p, 0, z) with the z-axis oriented along the magnetic field. We prescribe to each 
electron a definite value of the magnetic quantum number m^. Each single-electron 
wave function ^ depends on the variables and (p, z) 

^(p,0,z) = (27r)- 1 / 2 e-^M(^P) (2) 

where p denotes the numbering of the electrons. The resulting partial differential 
equations for ^(z , p) have been presented in ref. [20]. 

The one-particle equations for the wave functions ip^z, p) are solved by means of 
the fully numerical mesh method described in refs. [3,18,20]. In our first works on 
the helium atom in magnetic fields [18,20] we calculated the Coulomb and exchange 
integrals by means of a direct summation over the mesh nodes. But this direct 
method is very expensive with respect to the computer time and due to this reason 
we obtained in the following works [28-31] these potentials as solutions of the 
corresponding Poisson equation. The problem of the boundary conditions for the 
Poisson equation as well as the problem of simultaneously solving Poisson equation 
on the same meshes with Schrodinger-like equations for the wave functions ift^z, p) 
have been discussed in ref. [20]. 

The simultaneous solution of the Poisson equations for the Coulomb and ex- 
change potentials and Schrodinger-like equations for the wave functions ift^z, p) 
is a complicated computational problem, especially for atoms in strong magnetic 
fields. The problem consists in the different geometry of the spatial distribution 
of the electron density and the potentials correspondent to this density. In strong 
magnetic fields the distributions of the electronic densities are compressed towards 
the z axis and look like needles directed along the z axis. The equations for the 
wavefunctions can be solved in finite cylindrical domains as done in refs. [3,18,20]. 
For strong magnetic fields 7 >> 1 these domains can be rather small in the p 
direction. On the other hand, the potentials created by these charge distributions 
cannot have such a strongly anisotropic form and the Poisson equations for them 
must be solved on meshes with the distribution of nodes not very different for the 
z and p directions. This means some loss of the precision for the wavefunctions 
due to a decrease of the number of nodes in the area of a large electronic density. 
The most difficult problem, however, is the different asymptotic behaviour of the 
wavefunctions and potentials. The wavefunctions of the bound electrons decrease 
exponentially as r — > 00 (r is the distance from the origin). This simplifies the 
problem of the solution of the corresponding equations in the infinite space be- 
cause it is possible either to solve these equations in a finite domain Q (with simple 
boundary conditions V'lan = or 9ip/dn\ dn = 0) with negligible errors for domains 
of reasonable dimensions or otherwise to solve these equations in the infinite space 
on meshes with exponentially growing distances between nodes as r — > 00. The 
solutions of the Poisson equations for non-zero sums of charges decrease as 1/r as 
r — > 00. In result, every spatial restriction of the domain Q introduces a signif- 
icant error into the final solution. In some other mesh Hartree-Fock approaches 
developed for diatomic molecules (e.g. see [32,33]) this problem has been solved 



for finite Q by introducing special boundary conditions for the potentials obtained 
from the asymptotic behaviour of the potentials. This approach, being in princi- 
ple approximate, requires additional calculations with different extensions of Q to 
estimate the error. 

In the present approach we address the above problems by using special forms 
of non-uniform meshes [28]. Solutions to the Poisson equation on separate meshes 
contain some errors Sp associated with an inaccurate description of the potential 
far from the nucleus. However due to the special form of the function Sp(h) for 
these meshes (where h is a formal mesh step) the errors do not show up in the final 
results for the energy and other physical quantities, which we obtain by means of 
the Richardson extrapolation procedure (polynomial extrapolation to h — [3,41]). 
The main requirement for these meshes is not an exponential, but a polynomial 
increase of the mesh step h when r — > oo. Moreover, this behaviour can be only 
linear one, i.e. h^ 1 = 0(l/r) as r — > oo. The error of the mesh solution in this case 
has the form of a polynomial of the formal step of the mesh h — 1/N, where N is 
the number of nodes along one of the coordinates. In practical calculations these 
meshes are introduced by means of an orthogonal coordinate transformation from 
the physical coordinates x p to the mathematical ones x m made separately for p 
and z. And the numerical solution is, in fact, carried out on uniform meshes in the 
mathematical coordinates x m . The characteristic feature of these meshes consists of 
rapidly increasing coordinates of several outermost nodes when increasing the total 
number of nodes and decreasing the actual mesh step in the vicinity of the origin. 
Due to this property we call these meshes "Run away" meshes. To be concrete we 
present here two such meshes: 

1. The "Plain Poisson" mesh is generated by the coordinate transformation 

«p = Aj^r (») 

— oo < x p < +oo, — 1 < x m < +1, A is a constant. This simplest mesh of this group 
is near to the uniform ones near the origin (i.e. a plot of the distance between the 
neighbouring nodes contains a large horizontal section close to x p = 0) and then 
this mesh smoothly transforms to the "run away" behaviour for x p — > oo. 

2. The "Atomic Poisson" mesh 

x p = A ^A±3^L ( 4 ) 
1 x m 

(b > 0) allows obtaining more precise results for atoms at reasonable values of 
b < 1 due to a more dense distribution of nodes near the origin. In fact, this formula 
provides three different types of behaviour in three different domains: (a) A uniform 
mesh in a small vicinity of x p = 0. This behaviour provides absence of irregularities 
in the finite-difference representation of the Hamiltonian. (b) |x p | ~ ^l^ml " the 
quadratic expansion of the mesh, (c) h^ 1 = 0(l/r) as r — > oo. The distribution of 
nodes for not too big distances from the origin (b) given by the simple formula (4) 



is similar to well known "Lagrange meshes" for atoms and provide similar precision 
of results. (See e.g. [34] for a definition of the "Laguerre mesh" (a mesh with nodes 
at zeros of the Laguerre polynomials) as a "Lagrange mesh" suitable for systems 
with the Coulomb potential and [35] (22.16.8) for approximate formulas for the 
zeros of the Laguerre polynomials). 

The overall precision of our results depends, of course, on the number of mesh 
nodes and, if necessary, can be improved in calculations with denser meshes. The 
most dense meshes which we could use in the present calculations had 120 x 120 
nodes. In most cases Richardson's sequences of meshes with maximal number 
80 x 80 or 60 x 60 were sufficient. 



Ill THE STRUCTURE OF THE ATOMIC GROUND 
STATE CONFIGURATIONS FOR THE LIMIT 

7 — > oo 

In this section we provide some qualitative considerations on the problem of the 
ground states of multi-electron atoms in the high field limit. These considerations 
along with the well known electronic structure of the ground states at 7 = 
present a starting point for the combined qualitative and numerical considerations 
given in the following section. At very high field strengths the nuclear attraction 
energies and HF potentials (which determine the motion along the z axis) are small 
compared to the interaction energies with the magnetic field (which determines the 
motion perpendicular to the magnetic field and is responsible for the Landau zone 
structure of the spectrum). Thus in the limit (7 — > 00), all the one-electron wave 
functions of the ground state belong to the lowest Landau zones, i.e. m M < for all 
the electrons, and the system must be fully spin-polarised, i.e. s ZjX — — \. For the 
Coulomb central field the one electron levels form quasi ID Coulomb series with the 
binding energy Eb = -^1 for n z > 0, whereas -Eg (7 —> 00) - ► 00 for n z = 0, where 
n z is the number of nodal surfaces of the wave function crossing the z axis. In the 
limit 7 — > 00 the ground state wave function must be formed of the tightly bound 
single-electron functions with n z = 0. The one-particle binding energies of these 
functions decrease as \m\ increases and, thus, the electrons must occupy orbitals 
with increasing \m\ starting with m — 0. 

In the language of the Hartree-Fock approximation the ground state wave func- 
tion of an atom in the high-field limit is a fully spin-polarised set of single-electron 
orbitals with no nodal surfaces crossing the z axis and with non-positive magnetic 
quantum numbers decreasing from m = to m = —N + 1, where N is the number 
of electrons. In result, we have for the first 10 atoms and positive ions in the limit 
7 — > 00 the following structure of ground state configurations which is a simple 



substitute of the periodic law at very strong magnetic fields 
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We shall often refer in the following to these ground state configurations in the 
high-field limit as |Ojv) - The states |Ojv) possess the complete spin polarisation 
S z = —N/2. Decreasing the magnetic field strength, we can encounter a series of 
crossovers of the ground state configuration associated with transitions of one or 
several electrons from orbitals with the maximal values for \m\ to other orbitals 
with a different spatial geometry of the wave function but the same spin polar- 
isation. This means the first few crossovers can take place within the space of 
fully spin polarised configurations. We shall refer to these configurations by noting 
only the difference with respect to the state (0^). This notation can, of course, 
also be extended to non-fully spin polarised configurations. For instance the state 
ls 2 2p_ 1 3(i_24 i /'_35g_4 with S z = —2 of the carbon atom can be briefly referred to as 
I Is 2 ) , since the default is the occupation of the hydrogenic series Is, 2p_i, 3d-2, ■ ■ ■ 
and only deviations from it are recorded by this notation. 

IV GROUND STATE ELECTRONIC 
CONFIGURATIONS OF ATOMS AND POSITIVE 
IONS AT ARBITRARY FIELD STRENGTHS 

Currently the carbon atom is the most complicated system with a thoroughly 
investigated structure of its electronic configurations for arbitrary magnetic fields 
and we start this section with a consideration of this atom. 

In the case of decreasing the magnetic field strength from very large values to 
7 = the fully spin-polarised ground state configuration of a multi-electron atom 
must undergo one or several crossovers to become finally the zero-field ground 
state configuration. This configuration for the carbon atom corresponds to the 
spectroscopic term 3 P. In the framework of the non-relativistic consideration this 
term consists of nine states degenerate due to three possible ^-projections of the 
total spin S z = —1,0, 1 and three possible values of the total magnetic quantum 
number M = —1,0,1. For very weak magnetic fields it is reasonable to expect 
values S z = —1 and M = — 1 for the ground state which can be described in our 



notation as ls 2 2s 2 2p 2p_ 1 . 

Thus the possible ground state configurations of the carbon atom can be divided 
into three groups according to their total spin projection S z : the S z — — 1 group 
(low- field ground state configurations), the intermediate group S z = —2 and the 
S z = —3 group (the high-field ground state configurations). This grouping is 
required for the qualitative part of the following considerations which are based on 
the geometry of the spatial parts of the one electron wave functions. 

We start our consideration for 7 ^ with the high-field ground state and sub- 
sequently consider other possible candidates in question for the electronic ground 
state for S z = —3 (see Figure 1) with decreasing field strength. All the one electron 
wave functions of the high-field ground state ls2p_i?>d-2^f-z5gSh-<i possess no 
nodal surfaces crossing the z-axis and occupy the energetically lowest orbitals with 
magnetic quantum numbers ranging from m = down torn = —5. We shall refer 
to the number of the nodal surfaces crossing the z axis as n z . The 6/i_5 orbital pos- 
sesses the smallest binding energy of all orbitals constituting the high-field ground 
state. Its binding energy decreases rapidly with decreasing field strength. Thus, 
we can expect that the first crossover of ground state configurations happens due 
to a change of the 6/i_5 orbital into one possessing a higher binding energy at the 
corresponding lowered range of field strength. It is natural to suppose that the 
first transition while decreasing the magnetic field strength will involve a transition 
from an orbital possessing n z = to one for n z — 1. The energetically lowest 
available one particle state with n z = 1 is the 2p orbital. Another possible orbital 




FIGURE 1. The total energies (in atomic units) of the states of the carbon atom as functions of 
the magnetic field strength considered for the extraction of the ground state electronic configura- 
tions with S z = —3. Our results (solid lines) and data taken from ref.21 (broken lines). Energies 
and field strengths are given in atomic units. 



into which the 6/i_ 5 wave function could evolve is the 2s state. For the hydrogen 
atom or hydrogen-like ions in a magnetic field the 2p is stronger bound than the 
2s orbital. On the other hand, owing to the electron screening in multi-electron 
atoms in field-free space the 2s orbital tends to be more tightly bound than the 
2po orbital. Thus, two states i.e. the ls2po2p-i3eL24/_35<7_4 state as well as the 
Is2s2p_ 1 3(i_24/_35(y'_4 configuration are candidates for becoming the ground state 
in the S z = —3 set when we lower the field strength coming from the high field 
situation. 

Analogous arguments lead to the three following candidates for the ground 
state in case of the second crossover in the S z = —3 subset which takes place 
with decreasing field strength: ls2s2p 2p_i3d_ 2 4/_ 3 , ls2p 2p_ 1 3d- 1 3d- 2 ^f-3 and 
ls2s2p_ 1 3d_i3d_ 2 4/_ 3 . It is evident that the one particle energies for the 3d_i 
and 2p obey E 3 d_ 1 > E 2po for all values of 7 since they possess the same nodal 
structure with respect to the z-axis and only the 3oLi possesses an additional 
node in the plane perpendicular to the z-axis. For this reason the configuration 
ls2s2p_i3d_i3d_ 2 4/_ 3 can be excluded from our considerations of the ground state. 
This conclusion is fully confirmed by our calculations. 

Similar reasoning given in detail in ref. [30] can be repeated for two other (S z = 
—2 and S z = —1) subsets of states and leads to the results presented in table 1 and 
in Figure 2. 

Figure 3 allows us to add some more informations to the considerations of the 
previous section. This figure presents spatial distributions of the total electronic 
densities for the ground state configurations of the carbon atom. More precisely, 
it allows us to gain insights into the geometry of the distribution of the electron 




FIGURE 2. Energies of the ground state configurations as a function of the field strength. Ver- 
tical dotted lines divide regions belonging to different Hartree-Fock ground state configurations. 



TABLE 1. The Hartree-Fock ground state configurations of the carbon atom in 
external magnetic fields. The configurations presented in the table are the ground 
state configurations at jxnia < 7 < 7max- 



no. 


I vain 


Tin ax 


The ground state configuration 


M 


s z 


#(7min) 


f 





0.1862 


ls 2 2s 2 2p 2p-i 


-1 


-1 


-37.69096 


2 


0.1862 


0.4903 


Is 2 2s2p 2p-i2p + i 





-2 


-37.9334 


3 


0.4903 


4.207 


Is 2 2s2p 2p-i3d_ 2 


-3 


-2 


-38.3359 


4 


4.207 


7.920 


ls 2 2p 2p_ 1 3d_ 2 4/_ 3 


-6 


-2 


-41.7369 


5 


7.920 


12.216 


ls 2 2p_ 1 3d_ 2 4/_ 3 5 5 _ 4 


-10 


-2 


-43.6397 


G 


12.216 


18.664 


ls2p 2p-i3GL 2 4/-355-4 


-10 


-3 


-44.9341 


7 


18.664 


CO 


Is2p_i3d_24/-35g_46/i_5 


-15 


-3 


-50.9257 



density in space and in particular its dependence on the magnetic quantum number 
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FIGURE 3. Contour plots of the total electronic densities for the ground state of the carbon 
atom. For neighbouring lines the densities are different by a factor of e. The coordinates z, p as 
well as the corresponding field strengths are given in atomic units. 



TABLE 2. The Hartree-Fock ground state configurations of the lithium atom in 
external magnetic fields. The configurations, presented in the table are the ground 
state configurations at 7 m i n < 7 < 7 max . 
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and the total spin. Thereby we can understand the corresponding impact on the 
total energy of the atom. The first picture in this figure presents the distribution 
of the electron density of the ground state of the carbon atom at 7 = 0. The 
following pictures show the distributions of the electronic densities at values of the 
field strength which mark the boundaries of the regimes of field strengths belonging 
to the different ground state configurations. For the high-field ground state we 
present the distribution of the electronic density at the crossover field strength 
7 = 18.664 and for three additional values of 7 up to 7 = 1000. 

For each configuration the effect of the increasing field strength consists in com- 
pressing the electronic distribution towards the z axis. However most of the 
crossovers of ground state configurations involve the opposite effect which is due to 
the fact that they are associated with an increase of the total magnetic quantum 
number M = Z)u=i m n- 

For the lithium atom the analogous arguments and calculations [29] lead to the 
scheme presented in table 2. 

V GROUND STATE ELECTRONIC 
CONFIGURATIONS IN THE HIGH-FIELD REGIME 

Let us consider now the series of neutral atoms and positive ions with Z < 10 
in the high field domain which we define here as the one, where the ground state 
electronic configurations are fully spin polarised (Fully Spin Polarised (FSP) regime 
S z = —N/2). The FSP regime supplies an additional advantage for calculations 
performed in the Hartree-Fock approach, because our one-determinant wave func- 
tions are eigenf unctions of the total spin operator S 2 . Starting from the high-field 
limit we investigate the electronic structure and properties of the ground states 
with decreasing field strength until we reach the first crossover to a partially spin 
polarised (PSP) configuration with S z = —N/2 + 1. 

The approach to this investigation is very similar to that described in the previous 
section and the picture of the crossovers associated with the ground states of atoms 
and positive ions A + is presented in table 3. It should be noted that, for atoms 
with Z < 6 and ions with Z < 7, the state | Is 2 ) becomes the ground state while 
lowering the spin polarisation from the maximal absolute value S z = —N/2 to 
S z = —N/2 + 1. For heavier atoms and ions we remark that the state |ls 2 ) is not 



the energetically lowest one in the PSP subset at magnetic field strengths for which 
its energy becomes equal to the energy of the lowest FSP state. For these atoms 
and ions the state \ls 2 2p ) is energetically lower than |ls 2 ) at these field strengths. 
For atoms with Z > 7 and positive ions with Z > 8 the intersection points between 
the state |ls 2 2p ) and the energetically lowest state in the FSP subspace have to be 
calculated. In result, the spin-flip crossover occurs at higher fields than this would 
be in the case of | Is 2 ) being the lowest state in the PSP subspace. In particular, the 
spin-flip crossover for the neon atom is found to be slightly higher than the point of 
the crossover \2p ) — |2p 3cLi), and, therefore, this atom has in the framework of the 
Hartree-Fock approximation only two fully spin polarised configurations likewise 
other neutral atoms and positive ions with 6 < Z < 10. It should be noted that 
the situation with the neon atom can be regarded as a transient one due to closeness 
of the intersection \2po) — |2po3cLi) to the intersection |2po) — |ls 2 2po)- This means 
that we can expect the configuration |2p 3cLi) to be the global ground state for the 
sodium atom (Z — 11). In addition an investigation of the neon atom carried out 
on a more precise level than the Hartree-Fock method could also introduce some 
corrections to the picture described above for this atom. 

Summarising our results we remark that the atoms and positive ions with Z < 5 
have one FSP ground state configuration |0jv) whereas the atoms and ions with 
6 < Z < 10 possess two such configurations \0 N ) and \2p ). 

Possessing total energies for all the atoms with Z < 10 we can compare these 
results with the adiabatic calculations [6,7]. Both our results (E 2 b) and calculations 
[6,7] (-Eid) are carried out in the adiabatic approximation. The difference between 
Eid and E 2 t> consists only in the usage of the adiabatic approximation for obtaining 




FIGURE 4. Relative errors of the total energy in the adiabatic approximation depending on the 
charge of the nucleus. Eid - the total energies in the adiabatic approximation (B12 = 0.1, 0.5, 1, 5 
[6], B12 — 2.35 [7]), i?2D ~ our two-dimensional mesh results. B i2 = J5/10 12 G. 



TABLE 3. Total energies (a.u.) of the neutral atoms and ions A + at the 
crossover points of the ground state configurations. 



z 


7 Atomic state(s) 


-_E(Atomic) 


Ionic state(s) 


-E(A+) 


2 


0.711 


<M, 


Is 2 ) 


2.76940 


|0jv) 


2.32488 


3 


2.153 


<M, 


Is 2 ) 


7.64785 


|0jv> 


7.00057 




2.0718 


Is 2 } 




7.65600 


|0jv), |ls 2 ) 


6.94440 


4 


4.567 


On), 


Is 2 ) 


15.9166 


|0jv> 


15.07309 




4.501 


Is 2 ) 


15.91625 


|0jv>, |ls 2 ) 


15.01775 


5 


8.0251 


On), 


Is 2 ) 


28.18667 


|0jv) 


27.16436 




7.957 


Is 2 } 


28.17996 


|0jv), |ls 2 ) 


27.10004 


6 


18.664 


On), \2 Po ) 


50.9257 


\0n) 


49.50893 




14.536 


2po> 




47.23836 


\0n), \2 Po ) 


45.77150 




12.351 


2p ) 




45.07386 


\2 P0 ), \ls 2 ) 


43.72095 




12.216 


2p ), 


Is 2 ) 


44.9341 


\ls 2 ) 


43.70075 


7 


36.849 


On), \2 P0 ) 


84.4186 


\0n) 


82.58182 




30.509 


2po) 




79.34493 


\0n), \2po) 
\2po), |ls 2 > 
\ls 2 ) 


77.41246 




17.429 


2po) 




66.72786 


65.26170 




17.398 


2p ), 


ls 2 2p ) 


66.69306 


65.25362 


8 


64.720 


On), \2 P0 ) 


130.6806 


\0n) 


128.4054 




55.747 


2po) 




124.1125 


\0n), \2 Po ) 


121.69825 




23.985 


2po), 


ls 2 2p ) 


94.3773 


\2po) 


92.78308 




23.849 


ls 2 2p ) 


94.3336 


|2p ), |ls 2 2p ) 


92.62502 


9 


104.650 


On), \2po) 


191.8770 


|0at) 


189.1446 




92.624 


2po) 




183.6944 


\0n), \2 Po ) 


180.7819 




31.735 


2po>, 


ls 2 2 Po ) 


128.1605 


\2p ) 


126.4414 




31.612 


ls 2 2 Po ) 


128.1125 


\2po), ls 2 2p ) 


126.2897 


10 


159.138 


On), \2 Po ) 


270.220 


|0iv> 


267.0112 




143.604 


2po) 




260.2740 


|0at), \2 Po ) 


256.8459 




40.672 


2 Po ), \ls 2 2p ) 


168.4734 


\2p ) 


166.6327 




40.559 


ls 2 2p ) 


168.4217 


\2 Po ), \ls 2 2 Po ) 


166.4863 



energies E 1T) instead of exact solution of the Hartree-Fock equations for E 2 n- Thus, 
the comparison of E 1D and E 2 b allows us to evaluate the precision of the adiabatic 
approximation itself and obtain an idea of the degree of its applicability for multi- 
electron atoms for different field strengths and nuclear charges. All our values lie 
lower than the values of these adiabatic calculations. It is well known, that the 
precision of the adiabatic approximation decreases with decreasing field strength. 
The increase of the relative errors with decreasing field strength is clearly visible 
in the table. On the other hand, the relative errors of the adiabatic approximation 
possess the tendency to increase with growing Z, which is manifested by the scaling 
transformation E(Z, 7) = Z 2 E(1,^ /Z 2 ) (e.g. [36,20]) well known for hydrogen-like 
ions. The behaviour of the inner electrons is to some extent similar to the behaviour 
of the electrons in the corresponding hydrogen-like ions. Therefore their behaviour 
is to lowest order similar to the behaviour of the electron in the hydrogen atom at 
magnetic field strength 7/Z 2 i.e. this behaviour can be less accurately described by 
the adiabatic approximation at large Z values. The absolute values of the errors 
in the total energy associated with the adiabatic approximation are in many cases 
larger than the corresponding values of the ionisation energies. 

VI MESH APPROACH FOR SINGLE-ELECTRON 
ATOMIC AND MOLECULAR SYSTEMS IN 
STRONG ELECTRIC FIELDS 

In this section we present some features of our approach With respect to its 
application to systems in strong external electric fields and correspondingly some 
relevant physical results. In contrast to the situation discussed in the previous 
sections atoms and molecules in external uniform electric fields have no stationary 
states, because for every state there is a probability that one or several electrons 
leave the system. Thus, when switching on the external uniform electric field, all 
the stationary states turn into resonances. Using the complex form of the energy 
eigenvalues 

E = E - iT/2 

one may consider quasi-stationary states of quantum systems similarly to the sta- 
tionary ones. In this approach the real part of the energy E is the centre of the 
band corresponding to the quasi-stationary state and the imaginary part T/2 is the 
half-width of the band which determines the lifetime of the state. In this communi- 
cation we consider systems which can be described by two-dimensional one-electron 
Hamiltonians. These systems include the hydrogen atom and the molecular ion 
in an electric field [37-39] and the hydrogen atom in parallel electric and magnetic 
fields [39,40]. The only electron present in such a system can leave it under impact 
of the external electric field. From the mathematical point of view the problem 
consists in obtaining solutions of the single-particle Schrodinger equation for this 



electron with the correct asymptotic behaviour of the wavefunction as an outgoing 
wave. Currently we have three different possibilities for fixing this asymptotics 
realised in our computational program: 

1. Complex boundary condition method. This method is described in detail in ref. 
[38] . The method is based on the fact that the single-electron Schrodinger equation 
for a finite system can be solved with the arbitrary precision in a finite area both 
for stationary and for quasi-stationary eigenstates. The case of stationary states is 
considered in [41,3]. The approach for the quasi-stationary states will be discussed 
following [38]. Figure 5 presents the potential curve for the simplest Hamiltonian 
of the hydrogen atom in an electric field 

H =--(— + -— + — -—) ---Fz (5) 
2 y dp 2 p dp dz 2 p 2 I r 



F is the electric field strength multiplied by the charge of the electron. Analogously 
to [40,3] the calculations can be carried out in an area Q which is finite along the 
direction z. For this coordinate we used uniform meshes. The boundary of the 
area z = —L z _ for z < (F > 0) (Figure 5) is determined from the condition of 
small values of the wavefunction on the boundary and, therefore, small perturba- 
tions introduced by the corresponding boundary condition [3]. The values of the 
wavefunction on the opposite boundary of the area (z = L z+ ) cannot be excluded 
from the consideration. We consider non-stationary states of the system decaying 
into continuum of free particles. In this process an electron leaves the system in 
direction z — > +oo and, thus, an outgoing wave boundary condition is to be estab- 
lished on z = L z+ . The form of this boundary condition can be derived from the 
asymptotic behaviour of the wavefunction for z — » +oo and has the form 




FIGURE 5. The potential energy for the hydrogen atom V(p — 0, z) in the external uniform 
electric field. 



where k = ^(E+Fz)] 1 ^ 2 is the wavenumber. Solving the Schrodinger equation with 
the Hamiltonian (5) and the boundary condition (6) established on a reasonable 
distance L z+ from the origin of the system we can obtain the complex eigenvalues 
of the energy and the wavefunctions of the type presented in the Figure 6. 

This straightforward approach enables obtaining precise results both for atoms 
and molecules from weak to moderate strong fields (for instance for the ground 
state of the hydrogen atom up to F = 0.20 — 0.25 a.u.). 

2. Classical complex rotation of the coordinate z in the form z — > ze l@ . In this 
approach we have obtained precise results for atomic systems in strong fields from 
the lower bound of the over-barrier regime up to superstrong fields correspond- 
ing to regime |Re.E| << |Ihl£7| [39]. On the other hand, this method cannot be 
immediately applied to molecular systems in our direct mesh approach [42,39]. 

3. Exterior complex transformation of the coordinate z. In our numerical ap- 
proach we transform the real coordinate z into a curved path in the complex plane 
z. This transformation leaves intact the Hamiltonian in the internal part of the 
system, but supplies the complex rotation of z (and the possibility to use the zero 
asymptotic boundary conditions for the wavefunction) in the external part of the 
system. The transformation can be applied both for atoms and molecules and pro- 
vides precise results for fields from weak up to superstrong with some decrease of 
the numerical precision in the regime |Rei£| << jlmE] [39]. 

The numerical results obtained by all three methods coincide And are in agree- 
ment with numerous published data on the hydrogen atom in electric fields (see 
e.g. [43-46] and references in [38]). 

Some of our results for the hydrogen atom in parallel electric and magnetic fields 
[39] are shown in Table 4. 




FIGURE 6. Real and imaginary parts of the wavefunction of the HJ" molecule in a longitudinal 
electric field F (a.u.). 



TABLE 4. The ground state of the hydrogen atom in parallel electric and magnetic fields 
at F= 0.1, 1, 5. 



7 


JP _ 

r — 


n 1 




JP 


1 

1 


JP _ 





E 


r/2 




E 


r/2 


E Q 


r/2 





-0.5274183 


7.26904( 


-3) 


-0.6243366 


0.6468208 


-0.1350071 


3.083929 


0.01 


-0.532390 


7.2624(- 


3) 


-0.629329 


0.646812 


-0.140005 


3.083925 


0.1 


-0.574600 


6.6392(- 


3) 


-0.673584 


0.646053 


-0.184739 


3.083589 


1 


-0.8443098 


9.5923(- 


5) 


-1.0421379 


0.577291 


-0.6077008 


3.050207 


10 


-1.7498730 






-1.9579187 


0.1173924 


-2.375552 


1.955678 


100 


-3.790110 






-3.8219215 


4.9988(-5) 


-4.3806709 


0.4909467 


1000 


-7.66247 






-7.66807 




-7.82561 


1.04701(-2) 



The second system which we present in this section is the hydrogen molecular 
ion Hg" in a strong longitudinal electric field. Our approach allowed us to carry out 
the first correct consideration of this system [37] and, in particular, to obtain the 
potential curves of its ground state, presented in Figure 7 (left). The minima in 
these curves give the equilibrium internuclear distances, presented in this Figure 
(right) as a function of the electric field strength. One can see that at a critical 
value of the electric field about 0.065 a.u. the minimum disappears and, thus, above 
this critical value the hydrogen molecular ion cannot exist. This critical value of 
the maximal electric field F c = 0.065a.u. = 3.3V/A for the molecule H~2~ is in a good 
agreement with experimental results by [47]. According this work molecule may 
exist in laser beam fields with intensity less than 10 14 W/cm 2 which corresponds to 
3V/A, and does not exist in more intense fields. 



-0.4 n E 

: 2.3 




R F 

FIGURE 7. Left - rotential curves for the ground state of the Hj" molecule in longitudinal 
electric field F (a.u.). Right - The equilibrium internuclear distance of the molecule as a 
function of the applied electric field strength. 



TABLE 5. Equilibrium internuclear distances and corresponding en- 
ergies and half- widths of the energy of the ground state of H J molecule 
in longitudinal electric fields. 



F 


Ro 


Eo r/2 


F 


i?o 


E 


r/2 


0.00 


1.997 


-0.60264 


0.04 


2.062 


-0.60686 


1.23(-19) 


0.01 


2.001 


-0.60289 


0.05 


2.112 


-0.60943 


7.32(-15) 


0.02 


2.012 


-0.60366 


0.06 


2.198 


-0.61285 


1.49(-11) 


0.03 


2.031 


-0.60497 


0.065 


2.28 


-0.61501 


3.94(-10) 



VII CONCLUSIONS 

In this communication we have presented a 2D fully numerical mesh solution 
method in its various applications to atoms and simple diatomic molecules in strong 
external electric and magnetic fields. Specifically these are calculations of atoms 
with Z < 10 and their positive ions in strong magnetic fields and the comprehensive 
investigation of the electronic structure of the ground states of the Li and C atoms 
in arbitrary magnetic fields. For the carbon atom seven different electronic ground 
state configurations for different domains of the magnetic field strength have been 
found. The investigation of the series of atoms with Z < 10 in very strong magnetic 
fields enables us to evaluate the applicability of the adiabatic approximation and 
to show its decreasing precision for heavier atoms. 

The mathematical technique developed for solving Schrodinger equations for 
quasi-steady states allowed us to obtain a series of results for the hydrogen atom 
in parallel electric and magnetic fields and for the ri~2~ ion in strong electric fields. 

Thus, the method described above allows us to obtain a number of new physical 
results partially presented in this communication. These calculations are carried 
out in the Hartree-Fock approximation for multi-electron systems and are exact 
solutions of the Schrodinger equation for the single-electron case. As the following 
development of the method we plan to implement the configuration interaction 
approach in order to study correlation effects in multi-electron systems both in 
electric and magnetic fields. 
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